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3 

** 


1.  INTRODUCTION 


The  infrared  signature  of  ships  and  their  wakes  is  most  often  described  in  terms  of  contrast  to 
the  background  infrared  signature  of  the  sea.  Since  the  contrast  signature  of  ships  and  their  wakes 
involves  the  signature  of  the  rough  sea,  an  important  task  for  signature  modeling  is  to  calculate  the 
radiance  of  the  sea  surface  as  a  function  of  viewing  angle.  The  complication  in  the  modeling  arises 
because  of  the  background,  which  is  a  rough  sea  surface.  Wakes,  on  the  other  hand,  appear  as 
smoother  surfaces.  Since  the  infrared  wavelength  is  much  shorter  than  the  length  scale  of  distur¬ 
bances  on  the  water’s  surface,  the  rough  sea  may  be  thought  of  as  an  ensemble  of  small  facets,  each 
of  which  emits  and  reflects  infrared  radiation.  Facets  can  reflect  radiation  from  the  sky  and  from 
other  facets.  In  cases  where  the  length  (corresponding  to  the  projection  of  the  instantaneous  field  of 
view  (IFOV))  of  the  sensor  onto  the  ocean  surface  is  greater  than  the  length  scale  of  sea  disturbances, 
reflection  and  emission  properties  may  be  calculated  by  averaging  over  facets  weighted  by  the  proba¬ 
bility  density  function  for  facet  orientation. 

The  distribution  of  facet  orientations  can  have  a  significant  effect  on  the  infrared  radiation 
emanating  from  the  sea  surface.  If  the  sea  were  smooth  and  flat,  all  the  facets  would  be  horizontally 
oriented.  However  when  a  rough  sea  surface  is  viewed  near  grazing,  most  of  the  facets  seen  by  the 
sensor  are  tipped  away  from  the  horizontal  toward  the  sensor  [Goodell  1971,  Menat  1976].  This 
phenomenon  impacts  upon  both  the  amount  of  radiation  emitted  from  the  sea  and  the  amount  of  radia¬ 
tion  reflected  from  the  sky.  The  average  tilt  angle  causes  radiation  from  elevation  angles  relatively 
higher  in  the  sky  to  be  reflected  as  would  be  the  case  if  the  sea  where  not  rough.  Since  the  strength 
of  sky  radiation  decreases  with  increasing  elevation  angle,  the  effect  of  sea  roughness  is  to  decrease 
the  amount  of  radiation  reflected  from  the  sky.  At  the  same  time,  the  decreased  angle  of  incidence  to 
the  facet  increases  the  emissivity  of  the  facet  [Schwartz  and  Hon  1986],  which  increases  the  amount 
of  radiation  emitted  from  the  facet  At  angles  very  near  grazing,  the  reflected  radiation  decrease  is 
larger  than  the  emitted  radiation  increase.  For  other  angles,  the  situation  is  reversed. 


la 
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A  major  theoretical  problem  is  that  some  facets  are  hidden  from  view  by  others.  This  is  called 
the  hidden  surface  problem,  or  shadowing  problem.  Two  approaches  to  this  statistical  problem  are 
possible.  The  first  approach  is  to  calculate  the  bidirectional  reflectance  distribution  function  (BRDF) 
under  the  assumption  that  there  is  no  shadowing  and  then  to  introduce  some  heuristic  corrections  to 
take  shadowing  into  account.  This  approach  has  been  studied  by  Torrance  and  Sparrow  [1967]  for 
the  case  of  rough  material  surfaces  and  by  Maxwell  et  al.  [1973]  for  the  case  of  rough  water.  The 
second  approach  is  to  calculate  the  contribution  to  the  total  surface  area  seen  by  a  sensor  made  by 
facets  of  a  given  orientation.  This  approach  permits  a  simple  adjustment  to  one  of  the  effects  of  sha¬ 
dowing.  1  he  second  approach  is  presented  in  this  report.  Both  approaches  lead  to  the  same  result  in 
cases  where  there  is  no  shadowing.  We  show  that  the  simple  correction  allowed  by  the  second 
approach  leads  to  reasonable  results  for  viewing  angles  near  grazing  incidence. 
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Section  2  presents  the  statistics  needed  to  generate  the  facet  distribution,  and  a  simple  model  is 
shown  to  account  for  shadowing.  The  orientation  of  each  facet  affects  the  emissivity  and  reflectivity 
since  both  have  angular  dependence.  Section  3  shows  the  relevant  geometry  and  how  facet  reflec¬ 
tivity  and  emissivity  depend  on  facet  orientation.  By  using  these  results.  Section  4  describes  how  the 
apparent  radiance  of  the  rough  sea  compares  with  that  of  a  smooth  sea.  Finally,  Section  5  presents 
the  application  to  the  problem  of  ship  wakes.  This  section  describes  the  conditions  under  which  a 
ship’s  wake  may  appear  to  have  a  temperature  that  is  warmer  than  its  background. 

2.  STATISTICS  OF  FACET  PROJECTED  AREAS 

Here  we  calculate  the  waveslope  probability  density  function.  An  analytic  result  is  available  for 
the  case  where  the  wave  amplitudes  are  sufficiently  small  that  linear  superposition  of  amplitudes  is 
valid.  The  essential  elements  of  this  calculation  have  appeared  in  the  literature  [Beckmann  and  Spiz- 
zichino  1963].  (For  completeness,  the  derivation  is  presented  in  Appendix  A.)  The  central  result  can 
be  expressed  in  terms  of  a  probability  density  for  the  orientation  of  a  unit  vector  n  normal  to  the  sur¬ 
face  of  a  facet.  The  simplest  description  of  this  density  is  in  terms  of  a  system  of  spherical  polar 
coordinates  with  the  polar  axis  (the  z  axis)  in  the  zenith  direction  and  azimuth  measured  with  respect 
to  the  upwind  direction  (the  x  axis).  If  0  and  0  are  (respectively)  the  polar  and  the  azimuthal  angles 
of  n ,  the  probability  density  is: 

p(6,<t>)  =  (1/2-tt)  ( ab)~l  tan0  sec2  0  e-‘an,(M«w,  (1) 


with 

cos  20  sin20 

«2  b2  J* 

where  a2  is  the  variance  of  the  derivative  of  the  waveslope  in  the  upwind  direction  and  b2  is  the  vari¬ 
ance  in  the  crosswind  direction.  Equation  (1)  may  be  thought  of  as  the  angular  probability  density  for 
finding  the  facet  at  a  given  place  on  the  surface,  oriented  such  that  its  normal  has  the  angular  coordi¬ 
nates  0,0.  We  assume  there  is  no  correlation  between  upwind  and  crosswind  directions. 

The  basis  of  the  approach  developed  in  this  report  is  the  calculation  of  a  fraction  of  the  total 
projected  area  /(0,0)  seen  by  a  sensor.  This  is  the  fraction  contributed  by  facets  falling  in  an  infini¬ 
tesimal  angular  range  center  at  0,0.  Once  this  fraction  is  determined,  it  can  be  used  to  weight  the 
reflection  and  emission  contributions  of  the  facets  to  the  net  observed  radiance.  We  define  two  types 
of  projected  areas.  The  first  type  is  the  horizontally  projected  area,  which  is  obtained  by  projecting  a 
facet  or  group  of  facets  onto  the  horizontal  plane.  The  second  is  the  line-of-sight  (LOS)  projected 
area.  LOS  projected  area  is  obtained  by  projecting  a  facet  or  group  of  facets  onto  a  plane  normal  to  a 
vector  parallel  to  the  LOS  to  the  sensor.  The  probability  density  function  p(0,0)  is  related  to  the  hor¬ 
izontally  projected  area  in  the  following  way:  the  horizontally  projected  area  AAh  contributed  by 
facets  with  orientations  in  the  infinitesimal  angular  range  A0A0  centered  at  0,0  is  given  by: 

AAh  =  Ahp(6,<t>)A6A<t>,  (2) 


-4(0)  =  j 


where  Ah  is  the  total  horizontally  projected  area  in  the  region  under  consideration.  The  unprojected 
surface  area  A A  of  the  same  group  of  facets  is  AA^/cos  (0).  That  A A  is  greater  than  A Ah  is  a  reflec¬ 
tion  of  the  well-known  fact  that  the  total  surface  area  of  a  rough  surface  is  greater  than  the  horizontal 
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f  area  that  it  covers.  The  LOS  projected  area  of  this  same  group  of  facets  A Aws  is  given  by 

A A  cos  (A).  where  A  is  the  angle  between  the  normal  to  the  facet  and  the  LOS  from  the  facet  to  the 
sensor.  The  expression  for  A Aws  in  terms  of  the  probability  density  function  is  thus: 

A Alos  =  Ahp(6,(f>)  ^5-22.  A0A<».  (3) 

cos  (0) 

The  total  LOS  projected  area  can  be  obtained  by  integrating  the  right-hand  side  (RHS)  of  Eq.  (3)  with 
respect  to  (6,4>)  over  the  upper  hemisphere. 

Notice  that  while  0  is  constrained  to  lie  between  0  and  it  12  there  is  no  similar  constraint  on  A. 
Consequently  the  RHS  of  Eq.  (3)  can  be  negative  for  some  orientations  of  the  LOS  to  the  sensor. 
Equation  (3)  associates  negative  projected  areas  with  facets  that  are  oriented  such  that  their  normal 
vector  points  away  from  the  sensor.  When  the  RHS  of  Eq.  (3)  is  integrated,  the  negative  projected 
areas  contributed  by  facets  that  point  away  from  the  sensor  partiaiiy  cancel  the  positive  projected 
areas  of  other  facets.  The  correct  contribution  for  these  facets  is  0,  since  they  face  away  from  the 
sensor.  (The  model  considered  here  does  not  take  multiple  reflections  into  account.)  We  modify  Eq. 
(3)  by  multiolying  by  H(ir/2  -  A),  where  H(x)  is  1  when  x  >  0,  and  zero  otherwise: 

A Alos  =  AhP(d,<t>)  h (ir/2  -  Q).  (4) 

cos  (0) 

The  preceding  projected  area  discussion  makes  it  possible  to  give  an  expression  for  the  frac¬ 
tional  projected  area  f(6,<t>).  The  quantity  /(0,</>)A0A4>  is  defined  to  the  fraction  of  the  LOS  pro¬ 
jected  area  seen  by  the  sensor  that  is  contributed  by  facets  with  orientations  in  the  range  A0A<f>  cen¬ 
tered  at  d,4>.  It  is  given  by 


P(6,<t>)  H (tt/2  -  A) 

_ cos  (0) _ 

( P(0,4>)  »(*/2  -  WM* 

J  cos  (0) 


(5) 


That  is,  is  the  expression  of  Eq.  (3)  is  divided  by  the  whole  projected  area  seen  by  the  sensor. 

The  limits  of  the  integration  here  and  throughout  this  report  are  understood  to  be  0  to  2ir  for  the 
azimuthal  angles,  and  0  to  ir/2  for  the  polar  angles. 

An  alternative  method  for  arriving  at  Eq.  (5)  involves  consideration  of  a  shadowing  factor.  This 
is  the  probability  that  a  facet  of  given  orientation  is  not  hidden  from  view  of  the  sensor.  This  prob¬ 
lem  has  been  analyzed  in  the  context  of  radar  detection  [Smith  1967].  Ignoring  correlation  of 
waveheight  and  slope,  the  result  is  that  the  shadowing  factor  is  proportional  to  //( tr/2  -  A).  The 
proportionality  factor  is  independent  of  facet  orientation.  The  step  from  Eq.  (4)  to  Eq.  (5)  can  be 
viewed  as  the  replacement  of  LOS  projected  areas  with  the  product  of  LGS  projected  area  and  the 
appropriate  shadowing  factor. 

The  central  idea  introduced  here  is  to  use  the  expression  of  Eq.  (5)  to  weight  facet  contributions 
to  the  net  radiance.  For  comparison  with  the  results  of  the  BRDF  approach,  notice  that  the  integral 
in  the  denominator  of  Eq.  (5)  can  be  explicitly  carried  out  in  the  nonshadowing  case.  In  this  case, 
the  function  H  is  identically  equal  to  unity.  Elementary  trigonometry  may  be  used  to  express  cos  (A) 
in  terms  of  (0,<£)  and  (0,<£): 

cos  (A)  =  cos  (0)  cos  (0)  +  sin  (0)  sin  (0)  cos  (<t>  -  d>) ,  (6) 
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where  0  and  4>  are  the  polar  and  azimuthal  angles  of  a  unit  vector  r  that  is  parallel  to  the  LOS  to  the 
sensor.  (See  Fig.  1  for  the  coordinate  system  specification.)  The  integral  in  the  denominator  of  Eq. 
(5)  thus  involves  two  terms.  The  first  term  evaluates  immediately  to  cos  (0)  because  the  probability 
density  is  normalized.  The  second  term  involves  the  integral  of  the  product  of  cos  ( <t>  -  <£)  and 
p{6,4>).  The  probability  functions  considered  here  have  inversion  symmetry,  i.e.,  p(0,<t>)  = 
p(6,i r  —  <t>).  With  this  condition,  it  is  easy  to  see  that  the  second  integral  is  0.  However  it  can  be 
shown  that  this  integral  vanishes  under  more  general  conditions.  The  integral  in  the  denominator  of 
Eq.  (5)  thus  evaluates  to  cos  (0)  in  the  nonshadowing  case.  Thus  in  this  case: 


/(0,4>)  =  p(6,<t>) 


cos  (fl) 

cos  (0)  cos  (0) 


(7) 


Z  (ZENITH) 


X  (UPWIND) 


Fig.  1  —  Local  coordinate  system  for  a  facet,  r  is  the  unit  view  vector,  r'  is  the 
reflected  unit  vector  from  the  sky,  and  n  is  the  unit  normal  to  the  facet. 


3.  EMISSIVITY  AND  REFLECTIV  ITY 

The  facet  emissivity  €  and  reflectivity  p  are  both  functions  of  the  angle  Q.  The  Fresnel  formula 
for  unpolarized  radiation  may  be  used  to  calculate  p(fl)  if  the  complex  optical  constant  for  water  is 
specified.  The  appropriate  value  is  obtained  by  computing  an  average  weighted  by  the  response  func¬ 
tion  of  the  sensor.  Kirchhoff  s  law  in  the  form  e  =  1  -  p  may  be  used  to  obtain  the  angularly 
dependent  emissivity.  The  intrinsic  radiance  R,  is  given  by  Rt  =  e(U)RB,  where  RB  is  the  radiance 
of  a  blackbody  at  sea  temperature.  The  reflected  radiance  Rr  is  given  by  Rr  =  p(Q)Rs(6'),  where  Rs 
is  the  sky  radiance.  The  sky  radiance  depends  on  the  angle  0'  between  the  unit  vector  r  and  the  ze¬ 
nith  direction.  The  unit  vector  r  is  parallel  to  the  ray  generated  by  the  reflection  of  the  LOS  from  the 
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sensor  by  the  facet.  It  has  the  angular  coordinates  0',4>'.  A  vector  expression  can  be  given  for  r  in 
terms  of  it  and  r : 

r '  =  2  cos  (f2)n  -  r  .  (8) 

An  expression  for  O'  in  terms  of  0  and  6,<t>  is  obtained  from  this  equation  by  taking  the  inner  product 
of  both  sides  of  the  equation  with  the  unit  vector  pointing  in  the  zenith  direction. 

cos  (O')  =  2  cos  (fi)  cos  (6)  -  cos  (0).  (9) 

4.  ROUGH  SEA  RADIANCE 

The  apparent  radiance  of  the  rough  sea  Ra  is  the  expected  value  of  the  sums  of  /?,  and  Rr .  This 
expected  value  is  given  by: 

Ra(h  =  fj2  do  d<f>f(0, </>){[  1  -  p(Q)}RB  +  p(fl)Rs(0')\.  (10) 


Here  the  form  of  /  is  given  by  Eqs.  (5)  and  (1),  0  is  given  by  Eq.  (6),  and  O'  is  given  by  Eq.  (9). 
The  Fresnel  formula  is  to  be  used  for  p(fi),  and  the  function  Rs  must  be  specified.  An  atmospheric 
modeling  code  such  as  Lowtran  VI  can  be  used  to  create  a  tabulated  form  of  this  function.  Since  /  is 
a  normalized  weighting  function,  it  can  be  verified  that  Eq.  (10)  satisfies  a  consistency  relation. 
Specifically,  in  the  case  where  Rs  is  a  constant  equal  to  RB ,  one  finds  that  Ra  =  RB.  This  result  is 
required  by  the  second  law  of  thermodynamics. 

As  mentioned  in  Section  1 ,  the  results  obtained  here^  agree  with  those  obtained  by  using  the 
BRDF  approach  when  there  is  no  shadowing.  In  this  case,  /  is  given  by  Eq.  (7).  If  this  form  is  sub¬ 
stituted  into  Eq.  (10)  with  the  p(0,<t>)  given  by  Eq.  (1)  with  a2  =  b2  =  a2  (isotropic  case),  the 
weighting  terms  in  Eq.  (10)  become: 

tan  (0)  cos  (fl)  exp  (-tan^flyo2)  ^ 

cos  3(0)  cos(0) 

Apart  from  notational  differences,  this  weighting  factor  is  the  same  as  the  one  in  Lecompte  [1973]. 

To  see  the  typical  magnitude  of  the  effects  caused  by  a  rough  sea,  a  function  Rs  must  be  chosen. 
As  an  example,  we  have  exercised  the  Lowtran  program  with  summer  midlatitude  parameters  and 
have  obtained  the  long  wave  infrared  (IR)  (8  to  12  pm)  sky  radiance  form.  Fig.  2.  For  convenience, 
the  value  of  the  blackbody  radiance  of  the  sea  Rg  that  appears  in  Eq.  (10)  was  taken  to  be  equal  to 
the  value  of  Rs  at  the  horizon.  This  corresponds  to  selecting  temperatures  of  the  sea  and  of  the  air 
just  above  the  sea  to  be  equal.  These  Lowtran  results,  together  with  Eq.  (10)  (with  a2  =  b2  =  a2/ 2) 
were  used  to  calculate  z(O')  =  1  -  Ra(0')/RB  for  a  set  of  values  of  a.  For  this  calculation,  the  opti¬ 
cal  constants  appropriate  to  the  long  wave  IR  were  used  in  the  Fresnel  reflectivity  formula.  Since  an 
empirical  relationship  exists  between  a  and  the  wind  speed  given  by  Cox  [1954]: 

a2  =  0.003  +  0.00512  w,  (12) 


where  w  is  the  wind  speed  in  m/s,  this  set  of  results  is  effectively  parameterized  by  the  wind  speed. 
Figure  3  presents  the  calculated  function  z(0')  for  selected  values  of  a.  For  values  of  O'  near  75° 
(15°  away  from  grazing),  modest  values  of  the  wind  speed  produce  substantial  effects.  The  sea  can 
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Fig.  2  —  Sky  radiance  as  a  function  of  the  cosine  of  the  angle 
from  the  zenith 

appear  to  emit  up  to  10%  less  radiation  than  a  calm  sea  would  appear  to  emit.  Effects  of  this  magni¬ 
tude  can  have  a  substantial  impact  on  ship  detectability. 

5.  APPLICATION  TO  SHIP  WAKES 

When  a  ship  moves  through  the  water,  it  leaves  two  distinctly  different  wakes:  Kelvin  wakes 
and  the  turbulent  core  wake.  The  wake  region  of  interest  here  is  the  surface  of  the  turbulent  core 
wake.  One  geometric  feature  on  the  sea  surface  that  distinguishes  the  turbulent  wake  region  from  the 
background  is  that  the  surface  of  the  wake  appears  to  be  relatively  smooth.  Many  of  the  features  of 
the  background  that  make  it  rough,  such  as  capillary  waves,  have  been  greatly  reduced  in  the  wake 
region.  Only  those  waves  having  very  long  wavelengths  survive.  Since  the  width  of  the  wake  region 
is  on  the  order  of  the  length  of  the  transom  of  a  ship,  the  surface  of  the  turbulent  wake  region  appears 
as  a  smooth  rectangular  strip  in  a  roughened  sea. 

In  the  case  of  the  turbulent  core  wake,  recent  IR  measurements  show  that  the  wake  regions 
appear  to  be  either  cooler  or  warmer  than  the  background.  Daytime  measurements  show  that  the 
wake  has  an  apparent  temperature  cooler  than  the  apparent  temperature  of  the  surrounding  rough  sea 
[Peltzer  and  Garrett  1987],  Other  models  of  the  wake  region  have  been  used  to  compute  the  differ¬ 
ences  in  temperature  between  the  wake  region  and  the  background.  For  example,  in  Stewart  [1987], 
by  using  a  Navier-Stokes  approach  to  the  turbulent  core  region  that  includes  thermal  diffusion,  both 
cooler  and  warmer  wakes  are  predicted  compared  to  the  background.  Cooler  wakes  are  predicted 
during  the  day,  and  warmer  wakes  are  predicted  during  the  night.  Clear  sky  conditions  are  assumed 
in  both  cases.  Peltzer  and  Garrett  [1987]  and  Stewart  [1987]  conclude  that  there  is  an  intrinsic  tem¬ 
perature  difference  at  the  sea  surface  between  wakes  and  their  backgrounds.  Measurements  were 
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Fig.  3  —  Apparent  change  in  radiance  as  a  function  of  angle. 
Each  curve  corresponds  to  a  value  of  a. 


made  by  Peltzer  and  Garrett  [1987]  with  their  sensor  at  near  normal  incidence  to  the  wake  region. 
Near  normal  incidence  (the  emissivity  of  water)  is  close  to  unity,  and  the  reflectivity  is  small.  As  a 
result,  measurements  made  at  nadir  consider  only  self-emission  from  the  surface,  while  reflections 
from  the  sky  off  the  water  into  the  sensor  are  minimized.  Here  we  argue  that  in  the  absence  of  a 
temperature  difference  between  the  wake  region  and  its  background,  an  apparent  contrast  may  still  be 
observed  between  the  wake  and  its  background. 

By  using  the  statistics  of  the  facet  surfaces  developed  in  Sections  3  and  4,  we  compute  an 
expected  contrast  between  the  rough  surface  background  and  the  wake  region.  Since  the  surface  of 
the  wake  has  been  smoothed  in  the  turbulent  core  region,  we  model  the  surface  of  the  wake  as  a  flat 
water  surface.  That  is,  the  probability  density  of  the  slope  in  the  wake  region  is  a  delta  function  at 
the  origin.  Equation  (5)  gives  the  density  used  to  compute  the  expected  radiance  of  a  rough  sea. 
Expected  contrast  is  defined  to  be  the  expected  radiance  of  the  smooth  sea  minus  the  expected  radi¬ 
ance  of  a  rough  sea.  For  a  fixed-view  angle  from  the  zenith,  the  contrast  between  the  wake  and  its 
background  is  computed  as  a  function  of  the  background  roughness,  or  a,  where  a  is  the  slope  vari¬ 
ance.  The  slope  variance  varies  linearly  with  wind  speed  and  is  given  by  Eq.  (12).  Figure  4  plots 
the  result.  The  units  of  the  contrast  are  W/m2/sr.  The  sea  radiance  and  horizon  radiance  used  in 
generating  Fig.  4  are  equal  and  have  a  value  of  34.48  W/m2/sr.  Results  for  several  view  angles  (75°, 
80°,  82°,  84°,  86°)  are  shown.  For  small  values  of  a  and  angles  ^82°,  the  contrast  is  positive  and 
increasing.  This  means  that  the  radiance  of  the  wake  region  is  greater  than  the  radiance  of  the  rough 
sea.  As  a  is  increased,  a  maximum  positive  contrast  is  reached  before  the  contrast  decreases.  Lower 
angles  of  view  (such  as  75°  and  80°)  have  negative  contrast  values.  Therefore,  for  viewing  angles 
sufficiently  large  (i.e. ,  near  grazing),  the  apparent  temperature  of  the  surface  in  the  wake  region  is 
warmer  than  the  apparent  temperature  of  its  background.  However,  the  intrinsic  temperatures  of  the 
wake  and  background  regions  were  equal. 
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Fig.  4  —  Contrast  between  the  wake  and  its  background  as  a  function  of  a.  Each 
curve  corresponds  to  a  different  sensor  angle.  The  wake,  sky,  and  background 
radiance  values  are  assumed  to  be  equal  to  the  value  of  34.48  W/m2/sr. 


To  analyze  the  major  competing  effects  that  contribute  to  the  contrast,  consider  a  single  facet  of 
the  rough  sea  where  the  surface  is  flat  but  with  unit  normal  located  at  angles  (0,0).  As  in  the  previ¬ 
ous  sections,  r  denotes  the  unit  vector  along  the  LOS,  and  r'  is  the  unit  vector  defined  by  Eq.  (8). 
The  spherical  angles  of  r  and  r'  are  (respectively)  (0,0)  and  (0',0')-  The  reflectivity  of  water  is 
denoted  by  p(ft),  where  ft  is  the  angle  defined  by  n  ■  r  =  cos  (ft).  The  slope  of  the  surface  facet  is 
given  by  tan  0.  The  contrast  radiance  at  view  angle  (0,0)  is  the  radiance  of  horizontal  surface  minus 
the  radiance  of  sloped  facet  at  (0,0)  and  is  given  by 

C(0,0;0,0)  =  /?B[p(ft)  -  p(0)] 

+ /?,(0)p(0)  - /?,(0')p(ft).  (13) 


Equation  (13)  includes  contributions  from  both  intrinsic  and  reflected  radiance  values.  The  angles  ft 
and  0'  are  given  as  functions  of  (0,0)  and  (0,0)  by  Eqs.  (6)  and  (7).  To  highlight  the  effects  of  the 
reflectivity  of  water  and  the  sky  radiance,  add  and  subtract  the  term  Rs(6)p{Q)  to  Eq.  13  to  get 

C(0,0;0,0)  =  [p(0)  -  p(ft)]  (R,(0)  -  RB) 

+  p(Q)(R,(0)  -  R,(0')).  (14) 


We  now  analyze  the  signs  of  the  first  and  second  terms  in  Eq.  (14).  If  the  angle  of  view  is 
close  to  grazing,  most  of  the  facets  seen  by  the  sensor  will  be  tipped  up  toward  the  sensor.  There¬ 
fore,  we  assume  for  the  moment  that  0  =  0.  Under  this  assumption,  cos  ft  =  cos  (0  -  0)  and 
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cos  (O')  =  cos  (0  -  26),  It'  6  and  0  arc  in  the  interval  (0,  ir/2).  then  cos  A  >  cos  0  and 
cos  (O')  >  cos  (6).  (Since  the  slope  of  the  sea  surface  is  small,  this  inequality  is  reasonable  for  a  suf¬ 
ficient!}'  large  9.)  Since  p  is  monotone  increasing,  p(6)  >  p(12>.  Note  that  Rid)  <  RB  for  all 
6  e  (0.  7T/2),  since  the  sky  radiance  is  equal  to  the  sea  radiance  at  6  =  tt/2  and  is  monotone  increas¬ 
ing.  Therefore,  the  first  term  in  Eq.  (14)  is  negative.  The  second  term  is  positive  since  the  sky  radi¬ 
ance  is  monotone  increasing. 

For  a  single  facet,  there  are  two  competing  effects  when  determining  the  contrast.  One  effect  is 
caused  by  the  difference  between  the  reflectivity  of  flat  water  and  the  reflectivity  of  rough  water. 
The  second  effect  is  caused  by  the  difference  between  sky  radiance  reflected  from  flat  water  and  that 
of  rough  water.  With  respect  to  surface  roughness,  the  sign  of  the  contrast,  therefore,  depends  on  the 
relative  change  between  the  reflectivity  in  the  first  term  and  the  sky  radiance  in  the  second  term. 
Before  analyzing  this  point  in  detail,  note  that  when  <t>  &  4>,  the  angle  inequalities  hold  when  <p  -  <j> 
is  sufficiently  small.  Computations  show,  however,  that  the  signs  of  the  expected  values  of 
p(6)  -  p(A)  and  Rs(6)  -  Rs(0')  do  not  change  if  6  is  sufficiently  close  to  grazing.  (See  discussion 
below.) 

Let 


and 


P (0,fcO.<t>) 


p(0)  -  p(A) 

P(0) 


Ps  (0,0;0,<£) 


Rs(0)  -  Rs(0’) 
Rs(6) 


(15) 


(16) 


Equations  (15)  and  (16)  define  the  relative  change  in  reflectivity  and  sky  radiance  as  a  function  of 
facet  normal  at  a  fixed  view  angle  6,(f>. 


We  claim  that  then  0  is  near  grazing  and  <f>  =  4>,  both  Pp  and  Ps  are  increasing  functions  of 
surface  roughness  or  slope  when  the  facets  are  tipped  toward  the  sensor.  If  the  slope  is  increased,  A 
is  decreased.  This  implies  that  P p  is  increasing.  Likewise,  increasing  the  slope  decreases  O' .  There¬ 
fore,  Ps  is  increasing  as  a  function  of  slope.  When  the  restriction  of  requiring  equal  azimuthal  angles 
is  relaxed  and  the  expected  values  of  Ps  and  Pp  are  computed,  they  are  both  still  increasing  functions 
of  a. 


By  using  the  relative  change  expressions  given  by  Eqs.  (15)  and  (16)  in  Eq.  (14),  we  find  that: 


C(M;M)  =  W ^((LdfM)  (RsC0)  -  Rb) 


where 


+  w2(e.<t>-,o,<t>)Rs(e), 


Ws(0S<t>\0,4>)  =  p(0)Pp(l4K0,ct>), 
^2(M;0,<*>)  =  p(Q)PsCo,4>,6,<t>). 


(17) 


(18) 

(19) 
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Equation  (17)  explicitly  shows  how  the  contrast  depends  on  the  weights  Wx  and  W2,  which  are  pro¬ 
portional  to  P p  and  Ps,  respectively.  Multiplying  Eq.  (17)  by  the  density  of  facets  seen  by  the  sensor 
and  integrating  ( 6 ,0)  over  the  hemisphere  yields  the  expected  contrast  in  terms  of  the  expected  values 
of  W |  and  fV2  at  a  fixed  view  angle  (0,0).  Figure  5  shows  the  compute^  results  of  the  expected 
weights  for  several  angles  near  grazing  as  a  function  of  a.  As  before,  we  assume  that  the  sea  radi¬ 
ance  and  the  sky  radiance  at  the  horizon  are  equal  and  take  values  of  34.48  W/m2/sr.  The  solid  lines 
are  the  expected  values  of  IF,,  while  the  dashed  lines  are  those  of  W2.  Notice  that  in  both  cases,  W x 
and  W2  are  increasing  functions  of  a.  As  a  is  increased,  Wx  has  a  greater  rate  of  increase  than  does 
W2.  In  particular,  the  ratio  W2fWx  is  on  the  order  of  0. 1  when  a  is  0.05.  The  implication  is  that  the 
reflectivity  at  near  grazing  angles  of  0  changes  more  rapidly  than  sky  radiance  as  the  surface  is 
roughened. 


Since  the  first  term  in  Eq.  (17)  is  negative,  the  relative  change  in  the  sky  radiance  needs  to  be 
greater  than  the  change  in  reflectivity  for  the  second  (positive)  term  in  Eq.  (17)  to  become  dominant 
If  we  examine  view  angles  near  grazing,  we  see  that  the  ratio  W2/W t  needs  to  be  small  from  the  fol¬ 
lowing  argument.  The  expected  contrast  at  (0,0)  is  nonnegative  if  and  only  if 


<  W2(6,4>)>  -  <Wl(e&)> 


Rb  -  *,(g) 
Rs(e) 


>  o, 


(20) 


where  <  >  denotes  the  expected  value.  Let  0  =  rr/2  -  e,  where  e  is  small  and  positive.  It  can  be 
shown  that,  to  0[e2|, 
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Since  the  sky  radiance  changes  slowly  near  the  horizon,  the  derivative  of  Rs  is  small  (less  than  unity), 
while  the  value  of  Rs  near  the  horizon  is  an  order  of  magnitude  greater.  Therefore,  for  small  positive 
contrast  values,  <W2>  must  be  on  the  order  of  0.1  e  <W/!>.  That  is,  the  sky  radiance  changes 
more  slowly  than  the  reflectivity  at  near  grazing  angles. 

Note  that  when  6  is  not  near  grazing,  either  one  or  both  of  the  expected  weights  may  be  nega¬ 
tive,  as  a  function  of  a.  For  example,  in  Fig.  4,  the  contrast  is  negative  when  0  =  75°.  In  this  case 
<  W x>  is  positive  and  increasing,  while  <  W2>  is  negative. 

This  discussion  shows  that  a  nonzero  contrast  between  a  wake  and  its  background  can  be  driven 
by  pure  reflection  instead  of  an  intrinsic  temperature  difference  between  the  wake  and  its  background. 
We  now  consider  wake  temperatures  that  are  cooler  than  background,  which  means  that  wake  black- 
body  radiance  is  lower  than  background  blackbody  radiance.  The  radiance  in  the  wake  region  is 
decreased  by  an  amount  equal  to  the  difference  between  wake  and  background  regions  reported  in 
Peltzer  and  Garrett  [1987],  The  sky  radiance  at  the  horizon  and  rough  sea  radiance  were  assigned  the 
same  value  of  34.48,  as  before.  We  perform  the  same  computations  as  in  Fig.  4  and  compute  the 
contrast  between  the  wake  and  its  background.  Figures  6  and  7  show  the  results.  In  Fig.  6  the  radi¬ 
ance  in  the  wake  region  is  33.93,  and  in  Fig.  7  the  radiance  in  the  wake  region  is  32.85.  Notice  that 
in  both  Figs.  6  and  7  the  overall  contrast  is  reduced.  Angles  such  as  82°  no  longer  have  a  positive 
contrast  for  all  values  of  a.  (Compare  Figs.  4  and  6.)  Figure  7  shows  that  angles  greater  than  86° 
exhibit  positive  contrast  only  if  a  is  large  enough.  Nevertheless,  given  an  intrinsic  temperature  differ¬ 
ence  where  the  wake  is  cooler  than  its  background,  we  conclude  that  it  is  still  possible  to  have  an 
apparently  warmer  wake  than  its  background  when  viewed  near  grazing. 


Fig.  6  —  Contrast  between  the  wake  and  its  background  as  a  function  of  o.  Same  as 
in  Fig.  4,  but  the  wake  radiance  is  cooler  than  the  background  radiance.  The  sky  and 
background  radiance  values  are  equal  to  34.48  W/m2/sr.  The  wake  radiance  is  equal 
to  33.93  W/m2/sr. 
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Fig.  7  —  Contrast  between  the  wake  and  its  background  as  a  function  of  o. 
Same  as  Fig.  6,  but  the  value  of  the  wake  radiance  is  less  than  the  value  used 
in  Fig.  6.  Here  the  wake  radiance  is  32.28  W/m2/s r. 


6.  DISCUSSION 

We  have  presented  a  model  of  the  contrast  between  a  ship’s  wake  and  its  background.  Model¬ 
ing  of  the  background  entailed  the  use  of  a  probabilistic  model  of  a  rough  sea,  while  the  ship’s  wake 
was  modeled  as  a  smooth  surface.  We  conclude  that  in  the  absence  of  an  intrinsic  temperature  differ¬ 
ence  between  the  wake  and  its  background,  the  apparent  temperature  of  the  wake  could  still  appear  to 
be  warmer  than  its  background.  We  observe  that  a  rough  sea  can  emit  a  substantially  lower  amount 
of  radiation  in  the  infrared  when  compared  to  a  smooth  sea  surface.  The  sea,  being  a  poor  emitter  of 
radiation  at  near  grazing  angles,  reflects  cool  sky  radiation  from  higher  points  in  the  sky.  It  therefore 
looks  colder  than  its  blackbody  temperature.  When  considering  near  grazing  angles  of  view,  the 
smooth  wake  reflects  radiation  from  the  sky  at  a  lower  angle,  and  therefore,  the  reflected  sky  radia¬ 
tion  off  the  wake  is  warmer  than  reflected  sky  radiation  off  the  rough  sea  background.  The  reflec¬ 
tivity  (and  emissivity)  of  the  surface  also  affects  the  contrast.  The  details  of  how  the  sky  radiance 
and  reflectivity  of  vhe  surface  affect  the  contrast  were  analyzed  in  Section  5,  where  change  in  reflec¬ 
tivity  is  compared  to  the  change  in  sky  radiance  as  a  is  increased. 

Although  real  temperature  differences  between  wake  and  background  were  reported  in  Peltzer 
and  Garrett  [1987]  and  Stewart  [1987]  where  the  wake  was  cooler  than  its  background,  the  wake  can 
still  appear  to  be  warmer  than  its  background,  as  described  in  Section  5.  However,  the  roughness  of 
the  surface  (a)  must  be  increased  and/or  the  angle  of  view  must  be  closer  to  the  horizon.  Increasing 
the  surface  roughness  decreases  the  reflected  radiation  from  the  rough  sea;  increasing  the  view  angle 
increases  the  reflected  radiation  from  the  smooth  sea.  Clearly,  if  the  sea  is  sufficiently  cool,  there  are 
no  angles  for  which  the  contrast  is  positive. 

Other  authors  have  modeled  the  rough  surface  phenomenon  as  well  but  without  the  aid 
of  a  probability  model  of  the  sea  slopes  [Borrego  and  Machado  1985  and  Wilf  and  Manor  1984].  In 
some  extreme  cases,  other  authors  rely  on  one  emissivity  value  of  the  sea  surface  for  all  angles  [Singh 
and  Warren  1983]. 
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Another  popular  but  realistic  method  of  computing  the  apparent  radiance  of  a  rough  sea  is  to 
create  a  surface  that  has  the  same  amplitudes  of  its  Fourier  components  as  a  rough  sea.  The  method 
is  called  the  power  spectral  approach,  and  it  has  been  used  either  directly  or  indirectly  in  simulating 
images  of  a  rough  sea  [Menat  1976,  Schwartz  and  Hon  1986  (for  a  description  of  the  technique). 
Chapman  and  Irani  1981,  Ben-Yosef  et  al.  1983,  Borrego  and  Machado  1985,  Sidran  1981,  Wilf  and 
Manor  1984],  The  power  spectral  method  is  computationally  intensive  since  a  large  number  of  facets 
must  be  considered.  Furthermore,  many  surfaces  must  be  generated  to  compute  statistical  information 
about  the  sea  surface  in  question  for  each  wind  speed.  Analyzing  the  power  spectral  model  is  diffi¬ 
cult.  Our  probabilistic  approach,  however,  immediately  generates  the  statistics  of  the  rough  sea  sur¬ 
face  with  a  minimum  amount  of  computational  effort.  The  physics  is  contained  in  the  model  at  a 
glance,  thus  making  it  easy  to  analyze  and  verify  in  certain  cases. 

Finally,  our  slope  magnitude  model  of  the  sea  surface  is  a  two-dimensional  model  based  on 
linear  superposition  of  waves.  Real  sea  surfaces,  however,  have  nonlinear  waves,  as  pointed  out  in 
Huang  et  al.  [1983,1984],  In  particular,  Huang  et  al.  have  derived  the  probability  slope  distribution 
for  a  nonlinear  wave  in  one  dimension.  In  Huang  et  al.  [1984],  a  small  parameter  ok  is  introduced 
where  a  is  the  standard  deviation  of  the  surface  elevation  and  k  is  the  wave  number.  By  using  the 
assumption  that  the  parameter  is  small,  the  probability  slope  density  can  be  shown  to  be  normal  to 
order  (ok)2.  However,  small  changes  in  wave  slopes  can  possibly  lead  to  substantial  changes  in 
reflected  radiation.  Therefore,  nonlinear  waves  may,  in  fact,  affect  the  contrast  when  viewed  near 
grazing. 
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Appendix  A 

DERIVATION  OF  SLOPE  PROBABILITY  DENSITY 


We  begin  by  considering  the  function  z  =  h{x,y)  defined  to  be  the  water  height  above  mean 
sea  level  at  the  position  (jr,y)  in  the  horizontal  plane.  The  function  h  may  be  regarded  to  be  a  ran¬ 
dom  function,  that  is,  h(x,y )  is  to  be  regarded  as  a  random  variable  for  each  value  of  x  and  y. 
Corresponding  to  this  type  of  random  variable  is  a  probability  density  function.  Because  the  surface 
of  the  ocean  is  translationally  invariant  in  a  statistical  sense,  the  density  does  not  depend  on  x  and  y. 
Accordingly,  we  denote  it  by  p.  The  value  of  h(x,y)  is  the  sum  of  contributions  from  waves  of  a 
large  range  of  wavelengths.  To  the  extent  that  the  contributions  of  the  different  wavelengths  are  sta¬ 
tistically  independent,  the  central  limit  theorem  applies  and  p  is  of  normal  form.  This  assumption  of 
independent  contributions  is  an  approximation  that  breaks  down  for  large  values  of  h  as  the  result  of 
nonlinear  hydrodynamic  interactions.  However,  it  is  a  good  approximation  for  our  purposes  here. 
The  wave  slope  is  the  gradient  of  h .  The  statistics  of  spatial  derivatives  can  be  calculated  by  the  use 
of  finite  differences,  namely: 


jto  _  h(x',y)  h(x,y)  (A1) 

dx  x'y 

for  x'  sufficiently  close  to  x .  (A  similar  relation  holds  for  dh/dy.)  It  is  established  that  wavelengths 
shorter  than  a  certain  cutoff  contribute  little  to  the  slope  statistics.  Consequently,  Eq.  (Al)  may  be 
considered  valid  as  long  as  |  x'  —  x  \  is  less  than  the  cutoff  length.  If  the  anisotropy  introduced  by 
the  wind  direction  is  ignored,  h(x' ,y)  and  h{x  ,y)  are  identically  normally  distributed  random  vari¬ 
ables.  They  are  correlated  with  a  correlation  length  depending  on  |  x'  -  x  | .  Since  the  difference 
of  two  correlated  normally  distributed  random  variables  is  itself  a  normally  distributed  random  vari¬ 
able,  it  follows  that  dh/dx  and  dh/dy  are  normally  distributed.  The  wave  slope  j  is  defined  to  be  the 
magnitude  of  the  gradient  of  h,  i.e.. 


s 


+ 


(A2) 


Since  dh/dx  and  dh/dy  are  normal,  independent  with  zero  mean  and  equal  variance,  the  random  vari¬ 
able  s  has  a  Rayleigh  density,  given  by 


ps{s)  oc  se  s2/o2 ,  (A3) 

where  o2  is  the  slope  variance.  It  is  more  convenient  to  work  with  a  description  in  terms  of  spherical 
coordinates  on  a  unit  sphere.  We  introduce  the  angle  coordinates  6,<t>  to  describe  the  direction  of  the 
unit  normal  n  to  the  facet  face.  The  angle  0  is  the  angle  between  n  and  the  normal -to-the-sea  surface 
(the  z  direction),  while  <t>  is  the  azimuthal  angle  defined  as  0  when  n  is  in  the  z  ,x  plane.  Because 
s  =  tan  6,  it  follows  from  a  change  of  variables  that 

p„(0,4>)  =  (1  /27r)2cr2  tan  0  sec2  6  e_tan29/<r!.  (A4) 
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It  is  easy  to  see  how  the  nonisotropic  case  in  which  the  upwind  and  crosswind  variances  are 
unequal  generalizes  he  probability  slope  distribution.  Assume  that  the  x  and  y  directions  correspond 
to  the  upwind  and  c  osswind  directions,  respectively,  and  also  suppose  that  the  mean  slopes  are  zero 
in  both  x  and  y  directions.  A  normal  vector  to  the  surface  at  (x  ,>> )  is  given  by 

N(x,y)  =  (d$/dx,d$/dy,  1)  =  [j2  4-  1]1/2  n(x,y).  (A5) 

Decomposing  the  vector  N(x,y)  into  spherical  coordinates  (6,4>)  on  the  sphere  of  radius  [s2  +  l]l/2, 
it  follows  that  sin  6  =  s/[s 2  +  1] 1/2  and  tan  6  =  s . 

Reasoning  as  before,  the  slope  density  p(dh  /dx  ,dh  /dy)  is  normal  in  both  arguments  with  stan¬ 
dard  deviations  a ,  a  and  correlation  coefficient  r .  We  make  the  change  variables  by  letting 

3<i>/3x  =  s 2  +  1  sin  6  cos  <j>  =  tan  6  cos  <j> , 

d<b/dy  =  s2  +  1  sin  6  sin  <f>  =  tan  6  sin  <t>.  (A6) 


Since  the  Jacobian  of  the  transformation  is  sec  20  tan  6,  the  new  density  in  terms  of  (0,<£)  is  given  by 


where 


pn(6,<j))  = 


sec  26  tan  6 


2irax oy  V(i  -  r2)  ***  \_2(l  -  r2) 


tan  26 


£W>,  r) 


.  cos  2<t>  ,  sin  2 <f>  2 r  cos  <f>  sin  </> 

£-(,<P,r)  =  - z —  +  - - —  - , 


<Jx<Jy 


(A7) 


(A8) 


and  n  is  N  /||/V||.  Notice  that  if  a  =  a  and  r  =  0,  then  the  probability  density  is  given  in  terms  of  6 
only  and  reduces  to  Eq.  (A4). 
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